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We investigate the entanglement properties of finite spin rings, with noncolhnear Ising interaction 
between nearest neighbours. The orientations of the Ising axes are determined either by the spin 
position within the ring (model A) or by the direction of the bond (model B). In both cases, the 
considered spin Hamiltonians have a point group symmetry, rather than a translation invariance, as 
in spin rings with coUinear Ising interaction. The ground state of these models exhibit remarkable 
entanglement properties, resembling GHZ-like states in the absence of an applied magnetic field 
(model B). Besides, the application of an homogeneous magnetic field allows to modify qualitatively 
the character of the ground state entanglement, switching from multipartite to pairwise quantum 
correlations (both models A and B). 

PACS numbers: 03.67.Bg,75.10. Jm,75.50.Xx 



I. INTRODUCTION 

Spin rings represent prototypical low-dimensional sys- 
tems with highly entangled ground states [H, [1| • In par- 
ticular, antiferromagnetic and isotropic exchange inter- 
action induces pairwise entanglement between nearest- 
neighbouring spins; moreover, it maximizes the concur- 
rence within the set of translationally invariant states 
with vanishing magnetization Q. Quantum entangle- 
ment in anisotropic Heisenberg models has also been in- 
vestigated, partly in relation to the separability of the 
ground state for specific values of the applied magnetic 
field [i'-'To'l. Indeed, the magnetic field can be used as a 
control parameter in order to engineer the ground state 
and thermal entanglement of the system. The inter- 
play between the system anisotropies and the field of- 
fer even wider possibilities if one assumes that this needs 
not be homogeneous, but can rather be controlled locally 

I, [Hill. 

Previous analyses were mainly devoted to systems with 
translational invariance, where the anisotropies in the 
spin-spin couplings are independent on the site. How- 
ever, physical implementations of low-dimensional spin 
systems are typically characterized by point-group sym- 
metries, rather than traslational invariance |13| . This is 
the case, for examples, of nanomagnets [l3 |. that repre- 
sent a rich class of molecular spin clusters, with widely 
tunable geometries and physical parameters. A number 
of ring-shaped nanomagnets has been investigated in the 
last years [l3|; some of these possesses attractive fea- 
tures for the encoding and manipulation of quantum in- 
formation 

EMI. 

In this paper we consider rings formed 
by equivalent spins, where the anisotropies in the spin- 
spin couplings reflect the point-group symmetry of the 
molecule and the local environment of each spin. In par- 
ticular, we focus on spin models with noncollinear Ising 
interaction between nearest neighbours [H, [l^ . The re- 
sulting spin Hamiltonians don't fall into any of the com- 
monly considered cases and - unlike the standard Ising 
model - can exhibit highly entangled ground states, also 
in the absence of an applied magnetic field. 




FIG. 1: (color online) Twisted-spin representations of an 
hexagonal spin ring, (a) In the case of model A, the com- 
ponents of each spin Sfc in the Hamiltonian {H^ , Eq. [l| refer 
to local reference frames, with the Xk axis pointing in the ra- 
dial direction. The general reference frame is chosen such that 
fi = X, being r*, the position of sj,. (b) In the case of model 
B, the local reference frames refer to each of the couplings 
between nearest neighbouring spins; the x and y components 
are thus defined for each bond, with || (rfc+i — r^). There- 
fore, the components of each Sfc in the Hamiltonian {H^ , Eq. 
[SI refer to two local reference frames, one for the coupling 
with Sfc_i and one for that with Sfe+i. 

II. THE MODEL 

In order to provide an intuitive picture of the spin 
ring symmetry, we introduce the anisotropic Heisenberg 
model in a twisted-spin representation. In particular, 
we consider the case where the spins s.^ are located at the 
vertices of a regular polygon, and the directions of the co- 
ordinate axes are determined either by the position (r^) 
of each spin within the ring (model A) or by the direction 
of the bond between the exchange-coupled spins (model 
B). In the case of model A [Fig. [l{a)], the Xi compo- 
nent of spin Sj is along the radial direction (f ^ = , with 
the origin O corresponding to the center of the polygon) , 
whereas the yi axis is defined so as to form with Xi a 
left-handed reference frame in the polygon plane. In the 
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FIG. 2: (color online) Energy spectra of the noncoUinear Ising 
Hamiltonians (blue) and (red) with A'' spins s — 1/2 
and (, = X,Y . For both model A and model B, the spectrum 
is independent on whether J^^ = ^ Jyy or Jyy = ^ J^^ 
and on the sign of the coupling. While the ground states 
of are twofold degenerate, the ground doublets of 
present the following splittings S (not appreciable in the fig- 



ure): 3/ J ~ 10~ 
tively. 



,2.0 X 10"*, 10"^ for iV = 6,8, 10, respec- 



case of model B [Fig. [TJb)], instead, the x components 
of two neighbouring and coupled spins, and s^+i, are 
both along the side of the polygon [xj || (ri+i — rj)]. 



A. Spin model A 

In model A, the anisotropics in the Heisenberg cou- 
pling reflect the position of each spin within the ring: 



N 



(1) 



fe=l 



where N is the number of spins and sat+i = si. The 
primed spin components can be expressed in terms of 
the general reference frame: 



cos0fe Sk.x + sm0fe Sk, 



cos (pk Sk'. 



Sm (pk Sk.: 



(2) 



where x^ = (cos 0fc , sin </)fe ) and = (— sin 0^, cos 0fe) 
in the xy basis, with (j>k = 2(fc — l)7r/iV. The Hamilto- 
nian can be thus rewritten in terms of the general 
reference frame: 



N 

E 



^J^*^ + Z?"'^^ Sfe,a Sfe-(-l,/3, (3) 



where the tensors J and D account for the symmetric and 

Jf - J'k 



antisymmetric components of the interaction {J^^ = jf " 



FIG. 3: (color online) Concurrence (red symbols, right axis) 
and squared modulus (blue symbols, left axis) of the overlap 
between the ground state of and the state \^^{(j>)) (with 
^ = F,AF and x~ ^y^') see Eq. I15[) . For even spin num- 
ber (squares), the values of p and C are independent on the 
direction and sign of the coupling, but different angles (j) are 
considered in the x = (</> = 0) and x = X {(j) = n/2) cases. 
In the case of odd spin numbers (stars) the F cases have been 
considered. 



and Df = 
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J^y)sin(0fc+i 



'fc+i + 0fe), 



(4) 



while Dxx = Dyy = 0. 

In the case = Jyy = J^, the couplings of the 
above Hamiltonian become site-independent, and they 
reduce to an isotropic symmetric contribution (J^^ = 
— J^ cos{2tt/N) plus an antisymmetric one {D^^ = 
J^ sm{2'K/N)). In the following, we shall focus on the 
noncoUinear Ising models Hy and Hy, corresponding to 

tA _ n J- tA ^ tA respec- 



'•^yy ~ ^ ^ "^xx — J and J^^ — ^ Jyy — J 



tivcly. Both and Hy include, in the general reference 
frame representation, both (anisotropic) symmetric and 
antisymmetric terms, with components that depend on 
the spin index. 



B. Spin model B 

In model i3, the preferential directions in the 
anisotropic Heisenberg model are determined by the di- 
rection of each bond: 



N 



tB If III 



ttB \ ' / jB II III I 7^ 

— / , [•JxxSk,xSk+l,x "T '-'yySk,ySk+ 



1..) • (5) 



fc=i 



Here, the components of each spin are referred to two 
local reference frames: in the couplings with the spins 
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Sfe_i and Sfe+i, these are (xfc_i,yfc_i) and {xk,yk), re- 



spectively 
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- R 

= Sfc • Xfc = 
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sin (^fc Sfc,y, 
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= Sfc • Yfe = 


cos(/7fc Sfe^y - 


sini^fc Sk^x, 


s'" 


= Sfc • = 


COS(^fc Sk+l^x 


+ sinipk Sk+i,y, 




= Sfc • yf 


COSipk Sk+l,y 


- SiTKfik Sk+l,x, (6) 


being xf 


= (cos sin 


^k), yf = (- 


- sin ipk , cos (/Jfc) and 



Vfc = ""/S + {'2k — l)7r/Af. In the general reference frame, 
the above Hamihonian can be written as in Eq. [3l where: 



jyy 

'^k 

•^k 
Dk" 



J^yCOS^ipk + jf^sinVfc, 
{Jxx - Jyy) cos Vk Sinifik, 
0. 



(7) 



In the case = Jyy = , the couphngs of 
become site-independent, and only the symmetric and 
isotropic contribution is retained: Jjf ^ = J^^ — , with 
J^y — = 0. In the following, we shall focus on 

the noncollinear Ising models and iJy , defined as 
Jyy = ^ Jf, = and J^, = ^ J^y = J^, respec- 
tively. These cases correspond to a symmetric exchange 
(D^^ = 0), with the principal directions that vary from 
one spin pair to another. 

We finally note that the above models A and B can- 
not, in general, be rephrased one in terms of the other. 
In fact, model A can be rewritten in the twisted-spin 
representation B only by adding an antisymmetric con- 
tribution: 



J. 

J, 
J 
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XX 
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- j4 sm\n/N) + J^y cos^tt/N), 



J4cos^(ViV) 
0, 



(44 + ^d)cos(ViV)sin(V7V). 



(8) 



Analogous considerations apply to the model B in the 
twisted-spin representation A. In this case, the equations 
can be obtained from the above ones by swapping the A 
and B apices, and by changing the sign in the expression 
of the anti-symmetric exchange coefHcient. 



C. Symmetry properties of the A and B models 

Both the A and B models belong to the D„/i point- 
group symmetry with n corresponding to the spin 
number (n = N). In fact, one can show that 
and are invariant under A'^ different 62 rotations 
exp(— • hkir/h), with J = J2iLii^i + ^i)- rota- 
tions have axes parallel to the spin positions, = 
with k = 1, 2, . . . , iV/2; the other N/2 rotations have 
axes that coincide with the bisectors of the polygon sides, 
rifc = {rk + rk+i)/\rk+rk+i\, with k = N/2 + 1, . . . , N . 
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FIG. 4: (color online) Quantum entanglement in an hexag- 
onal ring of s = 1/2 spins, model Hx with J^^ — —1, in 
the presence of an external magnetic field b — 6x. (a) Pair- 
wise entanglement, quantified by the concurrence (C) between 
the spins s^ and s;, with |A; — i| ~ l(black squares), 2 (red, 
multiplied by a factor 5), 3 (blue, factor 25). (b) Block en- 
tanglement between subsystems iSi and 52, (consisting of Ni 
and N2 — N — Ni consecutive spins, respectively), quantified 
by Tr(pi), with pi the reduced density matrix of <Si. Inset: 
Residual tangle of the single spin. 



The Hamiltonians and are also invariant un- 
der reflection [at) about the polygon plane xy, and un- 
der the Cn rotation e'xjp{—i2'KJz/Nh) around the ver- 
tical axis z. We note that the collinear XY model, 

HxY = Jxx J2k=li 1,1/'); belongs to 

the point-group symmetry, if the x' and y' axes co- 
incide with symmetry axes of the polygon defined by the 
spin positions, i.e. if x' and y' are parallel to two of the 
rife. In the following we shall assume for simplicity that 
this is the case, and in particular that x' — x and y' = y. 

As suggested by the high degree of degeneracy of its 
energy spectrum (see below), the collinear Ising mod- 
els Hx and Hy are also invariant under a number of 
additional transformations, that apply to the orbital or 
spin degrees of freedom separately. In the case of the 
hexagon (Fig. [Ij, for example, the symmetry operations 
of H^=x,Y include all the elements of the Dg/i group, 
where each transformation is applied to the orbital part 
only: 6 different C2 rotations exp(— iL • nkir/H), with 
L — J2i=i ^rid the 6 rotations axes rife defined as 
above; reflection about the polygon plane xy; Cn rota- 
tion exjp{~i2TTLz/Nh) around the vertical axis z, whih 
can be thought as the analogue of translational invari- 
ance for a system with periodic boundary conditions. If 
we assume for simplicity that the electrons occupy spher- 
ically symmetric orbitals centered at the positions , the 
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transformations of Dg/i simply result in permutations of 
the spin indices. The rotation of an angle 7r/3 around 
the z axis, e.g., induces the following transformations: 
) — > (sfe+i.x, Sfc+i,y, Sfc+1,2). In general, 
one can show that the Hxy Hamiltonian of an A^-spin 
regular polygon is invariant under transformations be- 
longing to the D„;i group (with n = N). The models A 
and i?, instead, are only invariant under reflection of the 
orbital degrees of freedom about the xy plane. 

As far as spin transformations are concerned, the 
collinear XY model belongs to the T)2h group, with the 
C2 axes that coincide with the x' and y' axes. The C2 
rotations thus correspond to exp(— iS'o.Tr/fi), with S = 
J2i=i — ^'1 y'^ ^- The collinear Ising model Hx 

(Hy) is additionally invariant with respect to continuous 
spin rotations around the x' {y') axis. The Hamiltonians 
Ha and Hb, instead, are only invariant under reflection 
about the xy plane, corresponding to the transformation: 
isk,x,Sk,y,Sk,z) — > i-Sk x,-Sk,y,Sk,z)- We finally note 
that the Hamiltonians and H^^i (^, ^' = X, Y) are uni- 

tarily equivalent: = U^^iH^iU^^ . In the case ic' = x, 

for example, Uxx = ®k=i^'^^i~^^k,z4'k/K)- Therefore, 
one can associate the symmetry operation U^^i C of H^i 
to any symmetry operation of C of (and vice versa) . 
As shown below, the collinear and noncollinear models 
become significantly different (i.e. no longer unitarily 
equivalent) in the presence of an applied magnetic field. 



III. RESULTS 
A. Energy spectra and trial wavefunctions 

The results presented below are based on the direct 
diagonalization of the Hamiltonians and for 
rings of s — 1/2 spins. The energy spectra of the 
noncollinear Ising models (Fig. [5]) are independent on 
whether J^x 7^ = Jyy [H^ ^ H^) or Jyy ^ = Jxx 
[H^ — i?y), and coincide with those of the collinear 
Ising model {H^^x,y)- As already mentioned in the pre- 
vious section, all these models are in fact unitarily equiv- 
alent. The twofold degenerate ground states of H^, like 

any other eigenstate IvP^), can be derived from those of 
H^i by the unitary transformation U^^i : 

m^)) =W«'I«J> = K=iR'M+^ + k7:)] |a^), (9) 

where RkAv) = e:xp{-isk,z^), 4> = 4>itC), and k = 
0, 1. The expression of \a^) depends on whether the cou- 
pling has a ferromagnetic (x = oi' an antiferromag- 
netic (x = AF) character: 

10 - It.;, •■•t.;.), (10) 

where | f.) and | ],x) are the eigenstates of the single-spin 
projection along the x direction. Finally, one can easily 



verify that = for ^ = whereas (j) = ±7r/2 if (^, ^') 
coincides with {Y,X) or (A, F), respectively. 

The spectrum of is characterized by a lower de- 
gree of degeneration with respect to that of (Fig. H]) , 
reflecting the lower symmetry of the former Hamiltonian 
with respect to the latter one. In particular, the ground 
state doublet presents a splitting S, whose magnitude de- 
creases with the number of spins, as reported in the figure 
caption. For each spin ring, the energy spectrum is inde- 
pendent on whether — or — Hy and on the 
sign of the coupling . In fact, all these Hamiltonians 
are unitarily equivalent, being: 

H^U-\ withZ^ = (g)f^ie-""''-^/2, (11) 
H^=U{-H^)U-\ withZ^ = ®^iJe-*"^^'=-, (12) 

where the latter equation also implies that the energy 
spectrum is symmetric with respect to the origin. In 
the following we thus refer, without loss of generality, 
to the case of with > 0. In all the considered 
cases, it was found that the ground state could be 
expressed as a linear combinations of a limited number 
of symmetry-adapted states: 

N/2 

I*?) = [®liRi,z i^i)] E E ^ri'i'r), (13) 

k=0 {vfc} 

where the total spin projection along z of each component 
is fixed by k (M = N/2 - 2k). The components l*^"), 
whose coefficients C^* are determined numerically, are 
given by 

n 

where the 2k elements 1 < < A^ of the vector spec- 
ify which spins are fiipped with respect to the reference 
configuration |af ) = | f^f^ . . . f^f^) . In Eq. [HI differ- 
ent vectors correspond to components that cannot be 
transformed one into another by rotating the spin ring 
of an angle 0; around the z axis. For example, the com- 
ponents of and can be represented by all the 
states where the only two down spins are nearest neigh- 
bours or next nearest neighbours, respectively. In the 
case Jx < 0, the state |af) in Eq. [T3]is replaced by its 
antiferromagnetic counterpart \af^). Besides, additional 
relations are found between the coefficients of 1*1'^'°) and 
l^7V/2-fc)' that are the spin-flipped versions of one an- 
other, depending on the model, direction and character 
(ferromagnetic or antiferromagnetic) of the coupling (Ta- 
ble m . The use of the expression [T3] as a trial wavefunc- 
tion for the ground state allows to reduce drastically the 
dimension of the Hamiltonian to be diagonalized, for ex- 
ample from 256 to 12 for A = 8 spins s = 1/2, or from 
1024 to 15 for A^ = 10. 
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Jxx ^ Jxx ^ Jyy ^ Jyy ^ 


model A 
model B 


-11 1-1 
1-1-11 



TABLE I: Ratio between the coefficients CZ*' and C, 



'N 11- 



in the ground state of the noncollinear Ising models B. The 
components I'J"^'" ) and |'!E>]^//2-fc) ^1- [Sla-re the spin-flipped 
versions of one another. 



B. Ground state entanglement without field 

Unlike the standard Ising model, the noncollinear one 
(H^) presents remarkable entanglement properties in the 
absence of an external magnetic field. In fact, we find 
that the ground state of approximately corresponds 
to a symmetric combination of the two degenerate ground 
states of H^: 



l*«(0)> 



1 

71 



\i3fm + \f3im 



(15) 



where = 7r/2 and £, = F = AF) for negative 
(positive). As reported in Fig. |3] (filled blue squares), 
the squared modulus of such overlap, p — 
increases with the spin number N and approaches 1 al- 
ready for TV = 10. As a consequence, the ground state of 
the noncollinear Ising model B essentially corresponds to 
a linear superpositions of macroscopically different com- 
ponents, where the state of each spin in |/3q ((/>)) is orthog- 
onal to the state of the same spin in \/3f{(f>)). In analogy 
with the magnetization and Neel vectors, whose coherent 
tunneling is expected to take place in the ground state of 
molecular nanomagnets with ferromagnetic and antifer- 
romagnetic Heisenberg interaction, respectively [l^ . we 
introduce here the vectors np and uaf, defined as: 



1 ^ 



k ^k,z\(Pk 



k=l 



)}(16) 



UAF {(f>) = J2 i^k + 4>) Sfc Rkz (0fc + <^)} . (17) 
k=l 

These vectors have maximum modulus and opposite ori- 
entations in the case of the two macroscopically differ- 
ent components in Eq. [13 \{l3^i(l})\n^i(t))\l3^{<j)))\ = 1 
and (/?«(</>) |n^(^)|/3«(0)) = -(/3«(</>)|n4 (</.)|/3«(0)), with 
^ = F, AF. We can thus summarize the result reported 
in Fig. [3] by saying that the ground state of the non- 
collinear Ising models presents coherent tunneling of 
the vectors n^(7r/2). 

Similar results can be found in the case of odd spin 
numbers. Here, the overall Hilbert space can be di- 
vided into two uncoupled subspaces, including either 
the states with M = —N/2 + 2{k — 1), or those with 
M = +N/2 - 2{k - 1) (being fc = 1, 2, . . . , N/2). This 
results in a twofold degeneracy of all eigenvalues. Be- 
sides, the ferromagnetic and antiferromagnetic cases are 
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FIG. 5: (color online) Quantum entanglement in an octagonal 
ring of s = 1/2 spins, model Hx with J^^. — —1, in the 
presence of an external magnetic field b — fax. (a) Pairwise 
entanglement, quantified by the concurrence (C) between the 
spins Sfc and sj, with |fc — /| = 1 (black squares), 2 (red, 
multiplied by a factor 5). (b) Block entanglement between 
subsystems Si and 52, (consisting of A^i and = N — Ni 
consecutive spins, respectively), quantified by Tr(pi), with pi 
the reduced density matrix of Si. Inset: Residual tangle of 
the single spin. 



no longer equivalent, and an additional degeneracy of 
the ground state is induced by spin frustration in the 
latter case. In Fig. [3l we report the squared modulus 
of the overlap (^^1$^) for odd N (stars), where 1*^) is 
obtained by diagonalizing within each of the above 
mentioned subspaces. The trend as a function of N re- 
sembles that obtained for even spin numbers. The pos- 
sibility of approximating the ground state of with 
the state |^'^(0)) applies also to spins s > 1/2. In these 
cases, the single-spin states | fx) and | Ik) that enter 
the definition of \a^) (Eq. [101) correspond to \mx = +s) 
and \mx = — s), respectively. For example, from analo- 
gous calculations performed on s = 1 spins (not reported 
here) the overlap between the ground state of and 
|'I'^(0)) is larger than for rings of 1/2 spins with equal 
N. 

From the point of view of quantum correlations, the 
states 1^^) are equivalent to the Greenberger-Horne- 
Zeilinger (GHZ) states, characterized by a genuine mul- 
tipartite entanglement and by vanishing pairwise entan- 
glement. Pairwise entanglement between nearest neigh- 
bouring spins is however present in the ground state of 
H§. In fact, the concurrence between nearest neighbours 
has finite values (red squares in Fig. [3]) , that decrease for 
increasing N, and tend to zero as the ground state tends 
to the GHZ-like state |*^(0)). The concurrence between 
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FIG. 6: (color online) (color online) Entanglement between 
two subsystems of consecutive spins (<Si and 52) quantified by 
Tr(pf ), being pi the reduced density matrix of <Si. The trace is 
computed for the ground state of + , and is displayed 
as a function of the number of spins of <Si, A''i — N — N2, 
in the case N = 8 and for different values of the magnetic 
field. Inset: Residual tangle (blue curve) and concurrence C 
(plotted as 1 — C, red) as a function of the magnetic field b. 



pairs of spins that aren't nearest neighbours (not shown) 
is zero in all the considered cases. 



C. Magnetic field induced entanglement 

The magnetic field can be used as a control parameter 
in order to tune the quantum correlations within the ring. 
In the case of model A, the two degenerate and separa- 
ble ground states |/?^ (</»)) can be coupled by applying an 

homogeneous in-plane magnetic field = foX^^i Sk,x- 
As a result, the degeneracy is removed, and ) tends to 
a linear superposition of macroscopically distinct states. 
The overlap between the ground state and |\E'^(0)) (Eq. 
[T5)l as a function of the magnetic field is reported in Fig. 
H (iV = 6) and Fig. [5] (TV = 8). The in-plane magnetic 
field reduces the symmetry of the system and breaks the 
equivalence between the noncoUinear Ising models 
and their coUinear counterparts, and that between the 
ferromagnetic and antiferromagnetic cases. One of the 
consequences of the symmetry reduction is the removal 
of the degeneracy in the ground state. In fact, the energy 
splitting S = El — Eq between ground and first excited 
states, increases with b (not shown). The field however 
also mixes the subspace spanned by |/3g (</>)) and \l3i{(j))) 
with additional components, thus reducing the overlap 
between j^-^^) and |$*(0)) (see Table HH. All this is 
reflected in the entanglement properties of the ground 
state. For low values of the field, the concurrence be- 
tween all pairs of spins is negligible [Figs. SJa) and[SlJa)], 
whereas any blocks of A^i consecutive spins is entangled 



with the complementary block of N2 = N — Ni (panels 
(b)). Moreover, the value of Tr(/9^) is close to 1/2, in- 
dependently on the partition (A^i). Finally, the residual 
tangle, that quantifies the genuine multipartite entan- 
glement and is defined as 4det(pfe) — J2i^k ^fk (with pk 
the reduced density matrix of the fc-th spin C'ik the con- 
currence between spins i and k), is maximized (insets). 
All these features are consistent with a GHZ-like form of 
the ground state. For high values of the field, the ground 
state tends to a factorizable form, as shown by the simul- 
taneous suppression of pairwise, block and multipartite 
entanglement. In the intermediate region (b ~ 0.4 for 
N = 6 and b ~ 0.3 for N = 8), the abrupt reduction 
of the residual tangle and of the block entanglement is 
accompanied by peaks in the values of the concurrence, 
not only between nearest neighbours. 

We note in passing that the noncoUinear Ising H^, 
combined with an homogeneous magnetic field, results in 
an Hamiltonian and in ground state entanglement prop- 
erties that are equivalent to those of a collinear Ising in- 
teraction in the presence of an inhomogeneous magnetic 
field, with radial field orientation at each spin site (b^ || 
Tfc). While the latter geometry might produce remarkable 
entanglement properties in mesoscopic (pseudo)spin sys- 
tems [ll|, the former one seems much more suitable for 
producing analogous effects in nanometer-sized objects, 
such as molecular nanomagnets. 

In the case of model B, we consider a magnetic field 
applied along the z direction, giving rise to an addi- 
tional term in the Hamiltonian: — Sk,z- Such 
field preserves the equivalence between the ferromagnetic 
i^xx < 0) ^'^d antiferromagnetic (J^ > 0) models, as 
well as that between the = and Jyy — cases. We 
find that, also in the presence of the field, the ground 
state of the noncoUinear Ising model B is well approxi- 
mated by a linear superposition of two macroscopically 
different states: 



1 



(18) 



where 



|7«(e,0)) = [^liRiA^i + cj^ + kir)] \al{9)}. (19) 

Here ^ = F or ^ = AF, depending on whether the cou- 
pling has a ferromagnetic or an antiferromagnetic char- 
acter: 



(0)) 



I AF 

la* 



{<=ii?^[(-i)''e]}inn.-.t.>- (20) 



Unlike |/?^(0)) and IPficj))) (Eq. El, the states 17^(6*, </>)) 
and \jf{d,(j))) are not mutually orthogonal, unless 6 — 
7r/2, so that |af (0)) = \a^{9)) and hii0,<P)) = |/3«(0)). 
In general, |(7q(0, 0)17^(0, 0)) = |cos0|^; the normal- 
ization constant in Eq. [TS]is thus C = 2(1 -I- | cos0|^). 
As detailed in Table HIl for increasing values of the field. 



7 



6(b|lx) 


0.2 0.225 0.25 0.275 0.3 


model 
A 


P 


0.956 0.947 0.936 0.913 0.0595 
0.5 0.5 0.5 0.5 0.5 


6(b||z) 


0.1 0.2 0.3 0.4 


model 
B 


P 

Om/tt 


0.985 0.980 0.960 0.909 0.862 
0.250 0.229 0.206 0.173 0.136 



TABLE 11: Square modulus p of the overlap between the trial 
wavefunction (S, (p)) a-nd ground state of + (upper 
lines) or H'^ + ijf (lower lines), with iV = 8. The angle 
Om is the value of 9 that maximizes p. In the case of the 
models A and B, the magnetic field is oriented along the x 
and z directions, respectively. The above results refer to Hx, 
with = —1; as to the model B, no difference emerges 
between Hx + Hi^ and Hy , nor between = +1 and 

- -1 

the overlap between l^*^) and l^'f </>)) decreases, while 
the value 9m of the angle 9 that maximizes such over- 
lap decreases. As 9m passes from 7r/2 (like in the case 
6 = 0) to lower values, the entanglement properties of the 
ground state deviate from those of a GHZ state (Fig. . 
In particular, the entanglement between a block of A'^i 
consecutive spins and the remaining 7V2 = N — Ni ones, 
quantified by Trjjof}, becomes an increasing function of 
the Ni , whereas for & = it is practically independent on 
TVi, as for a GHZ state. The pairwise entanglement be- 
tween nearest neighbouring spins, quantified by the con- 
currence, increases (red curve, figure inset). The multi- 
partite entanglement, as quantified by the residual tan- 
gle, decreases for increasing b (blue curve). Altogether, 
the perpendicular magnetic field thus induces a transi- 
tion from a predominantly multipartite entangled ground 
state to one with large pairwise quantum correlations. 



IV. CONCLUSIONS 



In conclusion, we have investigated spin rings that 
are coupled by noncollinear Ising interactions, whose 
anisotropy reflects the point-group symmetry of the sys- 
tem. The ground states of these Hamiltonians exhibit re- 
markable entanglement properties. In particular, in the 
case where the preferential directions for each spin are de- 
termined by the direction of the spin-spin bond (model 
B), the system ground state is characterized by 

a large multipartite entanglement and by a low degree 
of pairwise entanglement. In fact, the overlap between 
1*^) and a GHZ-like state increases with the number of 
spins iV, and approaches 1 already for = 10. A ver- 
tical magnetic field can be used to substantially modify 
such picture, enhancing the pairwise entanglement be- 
tween nearest neighbouring spins at the expense of mul- 
tipartite entanglement. In the case where the preferential 
directions for each spin are determined solely by its po- 
sition within the ring (irrespective of the bond direction) 
- model A - the noncollinear Ising Hamiltonian is uni- 
tarily equivalent to the standard Ising model, and thus 
the degenerate ground state doublet is spanned by two 
factorizable states. However, the application of an mod- 
erate (with respect to J) in-plane field splits such dou- 
blet and induces multipartite entanglement in the ground 
state l^'o )• For increasing values of the field, under- 
goes a sharp transition towards a separable ferromagnetic 
ground state, accompanied a peak in the pairwise entan- 
glement, not only between nearest neighbouring spins. 
While these results have been obtained for s = 1/2 spins, 
analogous behaviours emerge from preliminary calcula- 
tions performed with higher spin values. 
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